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ABSTRACT 

While of order a million asteroids have been discovered, the number in rig- 
orously controlled samples that have precise orbits and rotation periods, as well 
as well-measured colors, is relatively small. In particular, less than a dozen 
main-belt asteroids with estimated diameters D < 3 km, have excellent rotation 
periods. We show how existing and soon-to-be-acquired microlensing data can 
yield a large asteroid sample with precise orbits and rotation periods, which will 
include roughly 6% of all asteroids with maximum brightness / < 18.1 and lying 
within 10° of the ecliptic. This sample will be dominated by small and very small 
asteroids, down to D ~ 1 km. We also show how asteroid astrometry could turn 
current narrow-angle OGLE proper motions of bulge stars into wide-angle proper 
motions. This would enable one to measure the proper-motion gradient across 
the Galactic bar. 

Subject headings: gravitational lensing: micro — minor planets, asteroids 



1. Introduction 

Modern surveys are discovering asteroids in prodigious numbers, soon to exceed a mil- 
liorfl The applications of such discoveries range from study of the detailed structure of 
Solar System resonances to defense from (or at least evacuation in the face of) asteroid 
collisions with Earth. However, despite the exponential increase in asteroid detection, the 
number of asteroids with we ll-measured ro t ation periods and colors (as well as as orbital ele- 



ments) remains quite small. IWarner et al.l (120091 ) have assembled a catalog of minor planets 



with at least some rotation-period information, which in its most recent (November 2012) 
edition has 5877 entries^] Among these are 1113 main-belt asteroids with excellent period 
determinations (quality flag '3'; excluding '3-'). 



1 http://www. minorplanetcenter.net/iau/MPCORB. html 
2 http: / /www. minorplanet.info/lightcurvedatabase. html 
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Figure [T] shows three different views of this sample. The lower and middle panels 
show differential and cumulative counts as a function of estimated /-band magnitude at 
opposition. For this purpose, we assumed that inner, middle, and outer main-belt asteroids 
are at distances of 2.3, 2.9, and 3.4 AU from the Sun. We also assumed (V — I) = 0.7 (similar 
to the Sun). That is, / = H + 51og(x(x — 1)) — 0.7, where H is the the cataloged absolute 
magnitude and xAU is the assumed heliocentric distance. There is a clear break at about 
/ = 11, a rounding commencing at about I = 13, and a sharp drop-off for / > 16. The 
upper panel shows that this drop-off corresponds to asteroid diameters of D ~ 3 km in the 
inner main belt and about D ~ 20 km in the outer main belt. Hence, while orbital elements 
are being measured for vast numbers of asteroids, the rotation periods of small main-belt 
asteroids are not being probed by current studies. 

Observations of Near Earth Asteroids s how that objects do wn to diameters of D ~ 200m 



have a minimum rotation period of 2.2 hrs (IPravec et al.ll2007l ). presumably because a faster 



spin exceeds the gravitational binding force. At present there are insufficient numbers of 
main belt asteroids with measur ed rotation period s and D < 3 km to determine whether 



they also obey this spin barrier. iDermawanl (12004 ) attempted to add ress this quest i on by 



measuring the rotation periods of small main-belt asteroids. However IWarner et al.l (120091 ) 



concluded that none of the rotation periods are secure, i.e., have a quality factor '3'. 

The fundamental reason that only a tiny fraction (~ 10~ 5 ) of the small asteroids being 
discovered have good rotation periods is that the discovery photometry is very sparse and 
followup photometry is rare. By contrast, modern bulge microlensing surveys have quite 
dense sampling, often 10-25 epochs per night. Moreover, these high-cadence zones cover 
contiguous areas that span 4° or more, meaning that typical asteroids that intersect these 
fields spend several weeks within them. 

One fairly simple idea then, would be just to identify known asteroids passing through 
these fields and measure their rotation periods (and colors) using archival microlensing data. 
However, microlensing surveys can have considerably greater return than this for asteroid 
science. We show here that these surveys can, by themselves, completely characterize the 
orbits of a large, well-defined subset of asteroids I < 18.1 that pass through their high- 
cadence fields. This means that it is possible to assemble a sample of asteroids that have 
not only well-defined orbits, rotation periods and colors, but also well-defined selection. By 
contrast, the present sample of small-to- moderate size asteroids with both rotation periods 
and orbits is completely heterogeneous. 

Thus, present and future microlensing surveys could probe completely new regions of 
asteroid parameter space. 
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2. Asteroids in Microlensing Surveys 

Our focus here is the possibility of extracting a large statistically well-defined sample 
of asteroids from existing and planned microlensing data. By "statistically well-defined" , we 
mean primarily a sample whose orbital elements and selection procedure are well enough un- 
derstood that the sample can be rigorously compared with a population model that specifies 
the distributions of orbital elements and (solar system) absolute magnitudes. The (microlens- 
ing) data are taken completely without reference to the possibility that asteroids will be or 
have been detected. Therefore, the key question is: how well can the orbital parameters be 
measured from existing data? Hence, we will focus mainly on the question of whether orbital 
elements can be accurately measured from microlensing data alone for a substantial number 
of asteroids. 



2.1. Framing the problem 

Microlensing surveys target the Galactic bulge, with the center of high-cadence obser- 
vations being near (RA,Dec) = (17:58:00,-29:20:00), i.e., 6° south of the Winter Solstice. 
The size of the high-cadence region varies according to survey. In this study we will adopt 
the parameters of the OGLE- IV survey and will comment in Section [6] on how the results 
should be adjusted for other surveys. OGLE-IV surveys an area of 11 contiguous square 
degrees at least once per hour, with 3/8 of this area being observed three times per hour. 
Typical trajectories crossing this zone would intersect it for about 4°. This high-cadence 
region is embedded in a larger low-cadence area|f|, with one or more observations per night 
(May-August) over a contiguous area of about 50 deg 2 . These observations are carried out 
whenever the bulge is visible, except three nights per month, typically chosen to be those 
when the Moon is passing through the bulge. Seeing of 1" or somewhat higher is achieved 
on more than half of all nights. 

Figure [2] shows the geometry of highly idealized asteroids in co-planar circular orbits. 
The angles bulge-Sun-Earth and bulge-Sun- asteroid are labeled # e and 9 ast , respectively. By 
the law of sines 

sin e = x sin 6» ast (1) 

where xAXJ is the radius of the asteroid orbit. The transverse velocity of the asteroid relative 
to the Earth is then 

^rel = (fast ~ %)-L = ^ast COS 6> as t - % COS e . (2) 



3 http:/ /ogle. astrouw.edu.pl/sky/ogle4-BLG/ 
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Setting v re i to zero in this equation and combining it with Equation ([I]) and Kepler's Third 
Law yields the orbital phases at which the asteroids are stationary (the transitions from 
prograde-to-retrograde and retrograde-to-prograde motion) 



Note that for x = (2, 3, 4), we have 6* eiStat = ±(49°, 56°, 61°). Thus, to a good approximation, 
this defines a zone of retrograde motion during ±60 days of the Summer Solstice (21 June). 
We now argue that to a good approximation, asteroids will be discovered if and only if they 
pass through the high-cadence fields during this "retrograde season" . 

There are a number of factors that vary over the course of the "retrograde season". 
We will begin by discussing observability and its affects on photometric and astrometric 
measurements. In Sections [3] and H] we will examine the impact of additional effects, such as 
asteroid phase angle and degree of trailing. First note that at the Summer Solstice (and for 
15 days on either side) the bulge fields can be observed at airmass < 2 for 10 hours, centered 
on their transit of the zenith. At 60 days from Summer Solstice, they can be observed for 
about 4.5 hours on one side of transit and 2 hours on the other. Observations near Summer 
Solstice therefore have a modest advantage of a longer observing window, which we will see in 
Section [3721 scales as t^ s , i.e., a factor (10/6. 5) 3 / 2 = 1.9. However, near their turning point 
the asteroids are moving more slowly and therefore spend more time in the microlensing 
fields. This effect scales as v~^ 2 . Of course, real asteroids are not on circular co-planar 
orbits. They could not be observed 6° from the ecliptic if they were. But the point is these 
two factors go in opposite directions. 

On the other hand, most of the asteroids that "enter" the high-cadence fields just after 
the "retrograde season" has ended are in fact re-entering and were already discovered during 
the retrograde season. Similarly, those that exit just before the start of the retrograde 
season, will enter the high-cadence fields at a later point, during the retrograde season. 
So there are few additional discoveries during these periods. And further from the edges 
of the retrograde season, the bulge observation window is rapidly contracting, making orbit 
characterization difficult. Of course, some very interesting asteroids will be discovered during 
this period, including asteroids interior to the Earth's orbit and other difficult-to-detect 
asteroids. However, from the point of view of characterizing the main features of the survey, 
these are a distraction. 

Hence, since the great majority of discoveries will take place during ±60 days of oppo- 
sition (Summer Solstice), and since survey geometry does not vary much over this period, 
we will initially carry out our calculations at Summer Solstice, when the asteroids are at 
opposition. 




(3) 
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Finally we note that while Earth is traveling from — 0ast,stat to +0ast,stat> the bulge- field 
observations cover an arc 2# astjStat of an asteroid orbit. At the same time, the asteroids move 
(in the same direction) through an arc 2^ e stat x -3 / 2 . Therefore, one season of observations 
covers a fraction of all asteroids in these orbits of 

ft —ft t~3/2 

r ^ast^tat (y ®,stat" 1 ' /.\ 

/season • 
71 

For x = (2,3,4), we have /season = (0.027,0.029,0.028), i.e., roughly constant. Since, the 
"cross section" of the high-cadence fields is about 4°, this means that a 10-year survey of this 
type would detect roughly (0.028 yr" 1 ) (10 yr)(4°/20°) ~ 6% of all asteroids passing within 
±10° of the ecliptic that satisfy the magnitude limit (still to be established). 

Virtually all asteroids discovered in this fashion would have (V — I) colors from OGLE 
data, and lightcurves spanning of order 40 days (i.e., not just in the high-cadence regions). 
Since asteroid rotation periods are typically less than one day, these observations would be 
sufficient to measure the periods, as well as amplitudes of variation, with good precision. 
Hence they can also be used to identify binaries and measure asteroid shapes. 



2.2. Asteroid Detection 



Microlensing surveys use difference imaging flAlard fe Luptonl Il998l ; IWozniak! 1200 0) , in 
which a reference image is convolved to the seeing of each target image, and then subtracted 
from it, which removes all non-varying (constant) sources. All that remains are "difference 
stars" at locations where the flux has changed. These include microlensing events and 
other stars whose brightness has changed, but it also includes asteroids, whose positions 
have changed. Microlensing fields are quite crowded, which means that the asteroid will 
frequently overlap with field stars. However, the great majority of these stars are below 
the sky background and so will subtract out in the difference image without even adding 
appreciable noise above the general background. Hence, the signature of an asteroid in these 
microlensing fields is the same as it would be for a series of images taken of a high-latitude 
field, i.e., a moving object. The only difference is that bulge fields are more crowded, so 
occasionally the signal is degraded by confusion with bright or variable stars. 

In principle, one could think about finding asteroids that are too faint to detect in 
individual images by stacking images along (hypothetical) asteroid orbits. However, as we 
will show below, it is actually impossible to properly characterize the orbits of asteroids 
that are undetectable in individual images (using microlensing data alone). Since "typical" 
asteroids are moving at retrograde velocities v ~ 12kms _1 at distances r ~ 2AU from 
Earth, they move only ~ 30" per hour. Hence, in fields with observations at least once per 
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hour, it is straightforward to identify such tracks on individual nights. Associating tracks 
between nights is also not difficult. As we will show in Section 15. 2\ from a single night of 
data one can predict the position on the next night to within 0(3"), which is extremely 
precise compared to the surface density of asteroids. 

3. Orbit Determination 

3.1. Orbit Parameterization 

Defining an orbit usually means measuring the six Kepler parameters, i.e., the five in- 
variants of motion plus the orbital phase. However, from the standpoint of understanding 
the orbital errors, it is more convenient to study the mathematically equivalent parameteri- 
zation of six Cartesian phase-space coordinates at a given instant of time, i.e., three positions 
and three velocities. Since the Earth's orbit is known extremely precisely, these six coordi- 
nates are equivalent to the following six: 2 instantaneous angular positions, 2 instantaneous 
angular velocities, the instantaneous distance, and the instantaneous "radial velocity", all 
measured from the center of the Earth. If these six are specified, one can transform to Kepler 
parameters, and vice versa. 

Hence, suppose that one conducts observations over 20 days while an asteroid passes 
through the high-cadence fields and then considers the ensemble of Kepler models that are 
consistent with these data. In principle, these models might span a large range of Kepler 
parameters, but all will agree extremely precisely on the angular positions as seen from the 
center of the Earth at the instant that the field transits on each night of observations. This 
is because these quantities are very nearly direct observables. That is, any model that is 
consistent with the data must reproduce the arc of mean nightly positions of the asteroid, 
and hence will yield the same angular positions and velocities at a given instant. 

These precise nightly determinations have two distinct implications. First they imply 
that the angular positions and velocities at some fiducial instant near the midpoint of obser- 
vations are extremely well determined (four parameters). Second, they enable a particularly 
simple estimate of the measurement precision of the remaining two parameters, i.e., the 
distance and radial velocity at this same instant. This is because, conceptually, fixing the 
geocentric angular position and velocity enables one to determine the diurnal parallax offset 
from each individual astrometric measurement, rather than carrying out a fit to the entire 
data set. The asteroid-distance measurement on a single night is then a simple integral over 
these measurements. The standard error of the mean of these measurements then gives the 
error in the distance at the above fiducial instant, while the error in the slope is the radial- 
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velocity error. Hence we must begin by deriving the distance error from a single night of 
data given that we know the angular position and velocity at transit. 



3.2. Single-night distance measurements 

The angular displacement of the asteroid due to diurnal motion of the observatory is 
Aa = Ap/r where p = (p n ,Pe) is the projected position of the observatory relative to the 
center of the Earth and r is the asteroid distance. While p changes in both the north and 
east directions during the night, the amplitude of the changes in the north component is 
generally small and will be ignored here in the interest of simplicity. Hence, the precision of 
the parallax measurement 7r = AU/r to an asteroid at distance r from Earth derived from 
n measurements of astrometric precision a is given by 

r y/nvai(p e ) VlOy V0.5y V30mas/ V ; 

where ( e = [var(p e )] 1/,2 /i2® is the standard deviation of observatory motion relative to the 
Earth center in the east direction normalized to Hence, for example, the distance of an 
asteroid at r = 2 AU would be measured to a precision of 8.8 x 10" 3 AU, assuming the above 
fiducial parameters. Note that for observations spanning t ohs and centered on transit 



, 1 - sinc(27rt ob s/day) 
C e ~ y cos 5 y 0.55 (6) 

where 5 is the latitude of the observatory and where we have assumed that the asteroid 
declination is similar. The final evaluation assumes t Q bs = 10 hr and 5 = —30°. 

Note that in the regime of interest, 1 — sinc(y) ~ y 2 /Q oc t 2 bs . Since n oc t bs, this implies 
(j(7r) oc t ob 3 / 2 as claimed above. More detailed calculations show that this relation remains 
approximately valid for a time series whose center is somewhat offset from transit, such as 
those at the limits of the retrograde season. 



3.3. Radial velocity measurement 

A typical asteroid at opposition moves at about 0.2° per day. Hence it remains in a 4° 
field for about T Q b s = 20 days. Assuming a fraction / g d of the nights have good seeing and 
low background, and that these m = / g d^obs/day nights are roughly uniformly distributed 
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over T bs, then the mean radial velocity will be measured with precision 



12 a(ny (n yl/2 aQ ..J ,-1/2. Tobs .-3/2, r , 

(7) 

By contrast, the error in the mean distance (over the measured trajectory) will be a(r) = 
m _1 / 2 o"(7r)r 2 /AU. Hence, the ratio of these two errors will be 

^KOA;© f— au 



/— AU / T obs n- 1 

v 7 !^ — — — = 9.9( — ¥—) . 8 

v^T nhK V20davs/ v ; 



<r(f)/AU v e T ohs ' V 20 days 

Thus, for typical asteroids at a ~ 3AU, a"(tv)A> as t ~ 17a(r)/a. This implies that the radial- 
velocity error is by far the dominant error in the problem. Relative to this error, the five 
other phase-space coordinates are known with essentially infinite precision. 



3.4. Impact on orbital period estimate 

In the above approximation that all instantaneous Cartesian parameters except v r are 
known perfectly, the specific potential energy is also known perfectly, while the fractional 
uncertainty in the specific kinetic energy is 2v r a(v r )/v 2 st , where now v r is taken to mean the 
asteroid velocity in the direction of the Earth but the frame of the Sun. Hence, the fractional 
error in the orbital period P is 

V{ P ) = 3 <t(o) = 3 <7fl.E|) ^ vMVr) = Q Q2g ojVr) Vr/Vast f ^ast \ ~ l (Q) 

P "2 a 2 \E\ ~ w 2 st ' ' O.Slkms- 1 0.2 UTkms^ 1 / ' {) 
Thus, for the above fiducial parameters and a P = 5 year orbit, cr(P) ~ 50 days. 



3.5. Astrometric precision and systematic errors 

The above estimates were scaled to astrometric measurements of 30 mas. This is a typ- 
ical error for good-seeing, low-background OGLE-IV stellar images at I ~ 18.4 (R. Poleski, 
private communication 2012). However, asteroid images are not stellar, but trailed. We 
evaluate the impact of trailing in the Appendix in terms of r] = e/2a, where e is the length 
of the trail and a is the Gaussian width of the point spread function (PSF). This impact 
is more severe for astrometric measurements in the direction of trailing than the orthogo- 
nal direction. Because the astrometry is limited by diurnal parallax measurements, which 
are essentially East- West and therefore aligned with the direction of asteroid motion (and so 
trailing), we evaluate the astrometric degradation factor G^ 2 (t]) in the trailing direction. For 
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our fiducial geometry, i.e., observations at opposition of an asteroid at x = 3 and reasonably 
good seeing a = 0.5" and 100 s OGLE exposures, we find rj — 1 and therefore Gx(rj) = 0.8. 
Hence our fiducial precision is achieved at / = 18.1. 

One must also worry about systematic errors, the most important of which is differential 
refraction. Asteroids are 1-3 mag bluer in (V — I) than the red giants that mainly define 
the astrometric reference frame in bulge fields, and the scale of the differential refraction 
is ~ 7 mas per mag per airmass. Differential refraction is extremely important because the 
distance measurements in Section I3T21 were based on measuring the astrometric deflection as 
a function of position relative to the Earth's center. The latter is very strongly correlated 
with airmass. Now, differential refraction can be calibrated from stars of similar color to the 
asteroids. However, (R. Poleski, private communication 2012) finds a systematics error floor 
of c sys ~ 2-3 mas, even after correction. 

Since these errors are systematic, one cannot necessarily rely on reducing them according 
to "square-root of N" , as was assumed in the equations for a{f) and a(v r ). However, the 
error floor imposed by systematic errors is very different in the two cases. The measurement 
of er(f) will not further improve if a < ^/nm<j sys , whereas for a(v r ) the improvement stops 
for Co < \/na Bys . Since the orbit determination is fundamentally limited by the latter, 
this implies that statistical precision is limited at o~ = V^-Csys ~ 10 mas. Since this is 
substantially below the fiducial value adopted in Equation (JSJ), the impact of systematic 
errors is likely to be small. 



3.6. Precision of rotation periods and amplitudes 

Just as the astrometric precision is degraded by trailing, so is the photometric precision. 
However, this effect has the same magnitude as the astrometric degradation in the orthogonal 
direction, G l y /2 (l) = 0.93. 

At / = 18.4 typical OGLE photometry errors are a(I) < 0.05 mag. therefore, we adopt 
a(I) < 0.04 mag at I — 18.1 including the effect of trailing. Small asteroids that would be 
observed near this limit have (full) amplitudes A of 0.1-0.2 mag. To evaluate the period and 
amplitude errors, we adopt a simple sinusoidal model: I(t) = Iq[1+ (A/2) cos(27rt/P rot + 0)]. 
For N, roughly uniformly sampled data points over time T, the amplitude error is then 
[a(A)}- 2 = (N/T) Jdt(dI/dA) 2 or 

A 0.04V0.15/ \10J V0.5/ V20days/ v ; 

which will therefore typically be quite good. 
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Similarly, the fractional error (or rather, strictly speaking, the minimum variance bound) 



in the rotation period is given by [cr(P rot )] 2 = (N/T) J dt(dl / dP TOt ) 2 or 



<?(Pr< 



0.002 



Pot cr(I) ( A 



rot 



1 day 0.04 V0. 15 



\10/ \0.5/ V20days/ 



-3/2 



111 



However, while the amplitude can be estimated directly from the scatter (assuming 
that the measurement errors are known), the rotation period must be found by testing all 
possible folds of the lightcurve, via e.g., a periodogram or Fourier transform. This can lead 
to multiple minima, and hence unless a single period is decisively favored by these tests, the 
true uncertainty can be much larger than the minimum variance bound. Thus, real data 
must be fitted before the t rue period completeness can be assess e d as a function of P ro t, A, 
and cr(I). See for example I Warner &: Harris j201lh : lHarris et all hold ). 



4. Variation of Sensitivity 

For simplicity, we have so far evaluated the sensitivity of microlensing surveys under the 
highly idealized assumption that detections are made at opposition and that asteroids are 
typically at distance r = 2AU (x — 3). We now successively relax these two assumptions. 

As the "retrograde season" progresses from opposition toward prograde motion, there are 
four distinct changes that affect the survey sensitivity: 1) the observation window per night 
(airmass < 2) shrinks, beginning about 15 days after opposition; 2) the asteroid retrograde 
motion slows, so that it remains in the high-cadence fields for more nights; 3) the same 
slowing proper motion decreases the trailing and its consequent error degradation; 4) the 
asteroid phase angle increases, thus making it fainter. In addition, the distance to the asteroid 
changes, but this effect is relatively minor. These effects each impact the astrometric and 
photometric precision differently. Therefore, we examine these two impacts separately. 

As discussed in Section [21 the astrometric precision of combined parallax measurements 
scales as (f rc iAobs) 3 ^ 2 - Since both factors decline during the approach to prograde motion, 
they tend to cancel. However, in the approximation of co-planer orbits, the first term actu- 
ally goes to zero, which is both unphysical for extended observations and not actually true 
for non-co-planar orbits. To suppress this singularity we adopt orbital inclinations of 10°, 
which is typical for bulge ecliptic latitudes. We then find that the ratio (i> re iAobs) monotoni- 
cally declines (improving parallax precision) from opposition to the prograde boundary. For 
example, it falls by a factor ~ 1.6 at 40 days from the Solstice, implying a factor 2 improve- 
ment in astrometric precision. There is also an improvement from decreased trailing, but 
this improvement is modest because the degradation at opposition is itself minor. 
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However, these improvements are countered by declining brightness at non-zero phase 
angle a. To calculate this, we assume a flux factor 

Z=(l-g) exp[-A 1 (tan(a/2)) Bl ] + #exp[-A 2 (tan(a/2)) B2 ] (12) 



where A 1 = 3.33, A 2 = 1.87, B x = 0.63, B 2 = 1.22, g = 0.15 flDvmockl 120070 . We then find 
that the net impact of all three factors on astrometric precision is complex. There is overall 
degradation from to 35 days after opposition, which peaks at 25% about halfway through. 
Then there is increasing improvement toward the boundary of the retrograde season. 

The photometric precision scales inversely as the square root of the exposure time, 
~ (f rc iAorb) 1//2 , i-e., as the cube-root of the astrometric impact. Curiously, as we show in 
the Appendix, the impact of trailing on photometric precision also scales as the cube root 
of its impact on astrometric precision in the range of interest. However, the non-zero phase 
impacts astrometry and photometry equally. Thus, the overall adverse impact is worse, 
rising toward 0.45 mag at 35 days and then declining toward zero at the season boundary. 

These results imply that the magnitude limit for defining orbits (from microlensing data 
alone) will vary as a function of phase, up or down by about 0.3 mag relative the special case 
of opposition, calculated above. Moreover, the rotation-period and variability errors will on 
average be somewhat larger than those calculated in Section 13.61 

Finally we consider how changing the asteroid orbit impacts these precisions. We find 
that for other orbital radii, the variation with time from opposition is qualitatively similar to 
the case examined above for x = 3, so that the main differences are found by comparing the 
cases at opposition. Here there are two effects, both due to the higher proper motion for closer 
asteroids. First, the time spent in the high-cadence zone is reduced and second, the image 
trailing is more severe. Together, these nominally lead to a factor 2 decrease in astrometric 
precision at x = 2 relative to x = 3. However, because the asteroid is closer, less astrometric 
precision is actually required to determine its orbit. The most critical measurement is the 
radial velocity v r , whose error scales as r 2 . Hence, the orbit determination substantially 
improves at x — 2. 

There is also a photometric effect, but this is much smaller, roughly a 0.25 mag degra- 
dation at x = 2. However, it is important to keep in mind that asteroids of fixed physical 
characteristics are 9 times brighter at x = 2 than x = 3. 

In summary, the sensitivity of microlensing survey data to asteroids does vary with 
orbital phase and orbital radius. However, the variations are relatively modest within the 
"retrograde season" for main belt asteroids. And the key point remains: this sensitivity, 
while variable, is rigorously calculable for asteroids with any definite set of properties that 
one might consider. 
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5. Asteroid Recovery 

5.1. Orbital-period timescales 

The error in the orbital period given by Equation would be small enough to char- 
acterize the asteroid orbit for many purposes. But it would also mean that the asteroid 
could easily be recovered if it passed through the bulge fields during microlensing season on 
a previous or subsequent orbit. That is, at any given instant, there would be an essentially 
one-dimensional, roughly 0.028 * 360° ~ 10° track defining the locus of possible positions for 
reappearance of the asteroid. While there might be other asteroids near this track, the ten- 
tative identification of one such image would fix the orbit with essentially infinite precision, 
which would immediately lead to secure predictions of its position in other images. 

What fraction of asteroids could be recovered in this way? Recall from Section 12.11 
that about 2.8% of all previously discovered asteroids that pass through high-cadence bulge 
fields will be detected each year. Ignoring for the moment that the asteroids have a range 
of inclinations, this means that an average of 2.8% will be recovered each year. Of course, 
for asteroids on 5 ± 0.5 year orbits, none will be recovered during the first four years, and 
14% during the fifth year. Now, if an asteroid has an integer-year orbital period, then it 
would return to exactly the same place in the high-cadence microlensing field, regardless of 
inclination. However, if it is highly inclined and does not have an integer-year period, it will 
likely miss the high-cadence fields. But here it is important to point out that in contrast to 
the observations required for high-precision orbit determination, the recovery observations 
can be quite sparse. The regions of such sparse observations are much larger than the high- 
cadence regions. Moreover, such sparse observations over roughly 50 deg 2 were carried out 
from 2001-2009 by OGLE-III. Hence, it is likely that tens of percent of asteroids could be 
recovered using a combination of OGLE-III and OGLE-IV data. 

It is quite possible that the great majority of asteroids whose orbits can be characterized 
by microlensing-survey data will already be kno wn and can simply be m atched to entries in 



the Minor Planet Center database. For example, iPolishook et al.l ( 120121 ) carried out a blind 
survey using the Palomar Transit Factory, which covered a similar area to the high-cadence 
OGLE fields with 10 observations per night, for four nights. Of the 30 asteroids that they 
discovered with 18 < R < 18.5 (whose faint limit corresponds to I ~ 18.1), all were already 
known. In the next bin 18.5 < R < 19, only 2 out of 50 were previously unknown. 
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5.2. Single-night predictions 

For completeness, we estimate how well the reappearance of an asteroid can be predicted 
given a single night of uniformly-spaced data, with n observations over t Q b s centered at 
transit, each with cq precision. For this purpose, we approximate the projected motion of 
the observatory (relative to the Earth center) as R & cos 5 sin (ut) ~ R & cos 5[ut — (ut) 3 /6], 
where u = 27rday~ 1 . Then, if one naively ignores parallax, one will make an error in the 
geocentric angular velocity A/i ~ uR® cos 5/r. This implies that if the distance r is estimated 
with precision cr(r), then the error in will be 

cr(r) 

<r(/i) ~ uR® cos 5 — — . (13) 
Then, after some algebra, one finds 



"W- / 1M800 (^ b ,)-^. (14) 



r 2 V n Rq cos S 

Together, these imply that the error in the predicted position will be 



obs 

independent of distance r. 



m = day x ,0.) = J™™***, = 2.si^y ,2 (^Y 2 ^, (15) 
V ' 3 VW V n cut 2 ^ \10J VlOhr/ 30 mas' V ; 



6. Other Surveys 

We have carried out our estimates using the characteristics of the OGLE-IV survey 
(supplemented by OGLE-III for asteroid recovery). There are two other current microlensing 
surveys, MOA and Wise, and one planned survey that is well on its way toward commission- 
ing, KMTNet. Because the MOA survey operates from New Zealand, it faces substantially 
more difficult observing conditions than the OGLE survey does in Chile. The seeing is typi- 
cally worse and the number of nights with good transparency during the whole night fewer. 
Combined, these two effects probably degrade the distance measurements by a factor of a 
few. Nevertheless, the current incarnation of MOA, with 2.2 deg 2 field and 1.8m telescope 
has been in operation since 2007. It covers 22 fields (almost 50 deg 2 ) every 50 minutes, which 
means that the asteroids remain in the MOA fields almost twice as long as in the OGLE 
fields. Hence many asteroids could be recovered. 



The Wise survey (jShvartzvald &: Maoa 120121 ) is carried out from Israel using a 1.0m 
telescope and 1 deg 2 camera. Because it is in the northern hemisphere, its observing window 
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is too short for effective nightly distance measurements. However, being in the northern 
hemisphere, Wise could make crucial observations that would complement a bulge asteroid 
survey carried out by either OGLE or MOA. Wise could observe a broad swath toward the 
Galactic anti-center at low cadence during the northern winter. To recover asteroids, it 
would only be necessary to make two 2-min observations on each of two successive nights. 
Hence a campaign of 3 months could cover about 1000 deg 2 , leading to recovery of a large 
fraction of asteroids over 5 years. 

The KMTNet survey will be carried out by three 1.6m telescopes (in Chile, South 
Africa, and Australia) each with 2° x 2° cameras. The survey plans to monitor four fields 
with a 10-minute cycle time. Hence, near the Summer Solstice it will observe each field 
roughly 120 times per night from the two good-seeing sites (Chile and South Africa). It 
could therefore plausibly reach 2.5 log(yT20/T0) = 1.4 mag deeper than the I = 18.1 limit 
estimated above. And, at fixed magnitude, KMTNet could measure rotation periods much 
better, if this proves to be a problem for lower-cadence observations. On the other hand 
fewer asteroids could be recovered from KMTNet data alone, because it does not plan to 
survey outside the high-cadence fields and because it lacks a long time baseline. 



7. Application to Bulge Science 



Precise astrometry of asteroids passing though the bulge would enable wide-angle as- 
trometry of the bulge at faint magnitudes. For example, it would permit determination of 
the proper-motion gradient across the Galactic bar, an important Galactic parameter whose 
measurement has never even been attempted. 



Sumi et all ( 120041 ) measured the proper motions of 5 million Galactic bulge stars with 
11 < I < 18 over 11 deg 2 using OGLE-II data, and R. Poleski (2013, private communication) 
is in the process of measuring proper motions over about 50 deg 2 using OGLE-III data. 
These catalogs enable a wide range of bulge science. However, they are intrinsically narrow- 
angle, i.e., tied to a local reference frame, and so do not permit wide-angle astrometric 
measurements. 

For example, if the bar is inclined to our line of sight by an angle a and is rotating as 
a solid body at frequency Q, then after some algebra one finds that the bar proper motion 
A/j,(£) at Galactic longitude £ (relative to £ = 0) is given by 

A/x(£) = (Q + /isgrA*) sin ^(cot a cos £ — sin £) — VLjj sin £(cos £ + cot a sin £) (16) 



where /is g rA* = —6.38 mas yr 1 (IReid &: Brunthalerl 120041 ) is the observed proper motion 
of SgrA* in the direction of Galactic rotation and VLu = U & /R ~ (10kms _1 )/(8kpc) = 
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0.26 mas yr -1 reflects the radial motion of the Sun U & toward the Galactic center. Hence, 
ignoring small terms, the gradient of the proper motion is dfi/di = + /isgrA*) cot a — 
It would be interesting to measure both this parameter combination and any deviations from 
this predicted simple behavior. 

There are several possible paths to transforming narrow-angle into wide-angle catalogs. 
One would b e simply to wait for GAIA. A second would be to tie the OGLE proper motions 



to Tycho-II (IHog et al.ll2000l ). This would require a special set of short exposures taken over 
several years because essentially all Tycho-II stars are saturated in standard OGLE data. The 
precision of this approach would be fundamentally limited by the relatively small number 
of Tycho-II stars (~ 200 deg -2 ), which would yield a precision of 0.20 mas yr _1 per deg 2 
patch (or 0.22 mas yr _1 if V < 10 stars were eliminated due to satur ation even in shortene d 



exposures). Third, one could tie local OGLE astrometry to UCAC-4 (jZacharias et al.ll2013l ). 
which in turn is tied to Tycho-II. It is not completely clear that UCAC reductions are 
effective in crowded bulge fields. Establishing this would require detailed tests. Moreover, 
this procedure would face the same fundamental limitations as the underlying Tycho-II 
catalog. 

Asteroids with well-defined orbits provide an alternate path to wide-angle proper mo- 
tions. For example, any asteroid found in OGLE-IV data and recovered in OGLE-III would 
have an extremely well measured period, as well as other orbital elements. Its predicted 
position relative to an ensemble of nearby clump giants from the previous orbit would there- 
fore be known to high accuracy for each observation, based on OGLE-IV observations plus 
the period determined from matching to one position of OGLE-III data. The displacement 
of the measured position relative to the same ensemble of clump giants in OGLE-III data 
would then give the relative angular displacement over the elapsed orbital period. 



Scaling from the iPolishook et al.l (120121 ) detection rate, and assuming 5 years of OGLE- 
IV searches, a 30% recovery rate from OGLE-III data, and roughly 25 observations for each 
within OGLE-III, this would imply roughly 11000 measurements with typical precision of 20 
mas, spread out over 50 deg 2 . Since the typical elapsed time is 5 years, this leads to proper 
motion precision of about 0.26 mas yr _1 per deg 2 , which is comparable to a calibration based 
on Tycho-II. 



8. Conclusions 



We have shown that existing microlensing surveys, particularly the OGLE survey con- 
tain a vast wealth of asteroid data from an almost completely untapped region of parameter 
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space. These data could be exploited to construct a catalog of asteroids with well-measured 
orbital parameters, precise rotation periods and amplitudes, and essentially perfectly known 
completeness down to / ~ 18.1. Currently, there are about a dozen such asteroids within a 
magnitude of this limit, and these do not have well-defined completeness properties. These 
same data can also be used to find asteroid binaries and find asteroid shapes. 
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Distinguished University Fellowship from The Ohio State University. We thank Scott Gaudi 
for insightful discussions. We thank David Polishook, Radek Poleski, and David Bennett for 
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A. Astrometry of Trailed Images 



For a general point spread function (PSF), F(x,y) that is normalized to unity, and 
assuming ph oton-limited noise, the astrometric error <7n in the n directions is given by (e.g., 



Gouldlll995f ) 



(NAFijAxAy) 2 
a 2 (NFij + B)AxAy 



(Al) 



where N is the total number of source photons, (Ax, Ay) is the pixel size, B is the background 
photon surface density, F^ is the mean of F in pixel and AFij is the change in this 

number under the impact of a displacement an. If the pixel scale is small compared to the 
structure of the PSF, then this can be evaluated using, 

(n • VF) 2 



a A = N dx dy- 



(A2) 



F(x, y) + B/N 

We now consider a Gaussian PSF of width er, that has been trailed by total length e. Since 
the pixel size is small, we can without loss of generality assume that the trail is in the x 
direction: 

ex P (-y 2 /2a 2 ) f^ 2 /2 d£exp(-(x + if /2a 2 ) 



F(x,y) 
from which we derive 



dF_ _ exp(-y 2 /2a 2 ) 
dx 2na i e 



2na 2 e 



-(x+e/2) 2 /2a 2 



-(x-e/2) 2 /2a 2 



and 



(A3) 

(A4) 
(A5) 



We now restrict to the background-limited case, i.e., F(x,y) + B/N — > B/N in the 
denominator of Equation fl A2|) . We then obtain for the trail direction, 



-2 



N 2 



1 — exp(— rf) 



e 



(A6) 



8tkj 4 B 



rf 




We recognize the first factor as the standard expression for the astrometric error of a Gaussian 
PSF below sky. The second term quantifies the degradation due to trailing. In the limit 
of T] <C 1, G x (rj) — > 1 — i] 2 /2, implying that the error is fractionally degraded by rf /A. For 
the intermediate case, rj — 1 (e.g., e = 1" and FWHM ~ 1.2") the error degradation is 
Gl /2 -> 0.80. For long trails, Gl' 2 -> rf 1 . 

For the non-trailed direction, the astrometric precision is 



Finally, we note that it is straightforward to show that the photometry precision of a 

1/2 

trailed image is degraded by exactly the same factor G y as the non-trail direction astrom- 
etry. 
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Fig. 1. — Lower and middle panels: Differential and cumulative number counts of asteroids 
with excellent rotation-period determinations as a f unction of estimated /-band magnitude at 
opposition as derived from the Nov 2012 update of I Warner et al.l (120091 ) . Apparent breaks in 
the distribution are indicated by dashed vertical lines: below / = 11, the distribution obeys 
a diameter distribution N oc D~ 2 (dot-dashed line), it begins to turn over at / = 13, and 
suffers a rapid drop-off for / > 16. This is compared to our estimated completeness limit 
from an OGLE-based survey of I = 18.1. Upper panel: estimated diameters of these same 
cataloged asteroids. Very few are smaller than 3 km. 
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Fig. 2. — Idealized geometry of asteroids passing through bulge micro lensing fields. As the 
Earth moves (orbital phase in its approximately circular orbit, the observations target 
the same bulge fields. The angle bulge-Sun-asteroid to an asteroid at xAU can therefore be 
calculated from the law of sines: sin 6^ = x~ l sin 0®. If the asteroid orbit is approximated as 
circular, then at 6 1 ® = 0®, s t a t = sin _1 [x(l + x + x 2 )~ 1 ^ 2 } the asteroids transition from prograde 
to retrograde motion (or vice versa). See Equation (j3J). Specifically, 49° < 0®, s tat < 61° for 
2 < x < 4. The season when many new asteroids can be discovered is approximately 20® )Stat . 
The fraction of asteroids discovered during this time is /season = [#ast,stat — #©,stat£~ 3//2 ]/7r, 
i.e., 2.7%-2.9% per season. 



